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The interaction between a single hole and a two-dimensional, paramagnetic, homogeneous electron 
gas is studied using diffusion quantum Monte Carlo simulations. Electron-hole relaxation energies, 
pair-correlation functions, and electron-hole center-of-mass momentum densities are reported for a 
range of electron-hole mass ratios and electron densities. We find numerical evidence of a crossover 
from a collective excitonic state to a trion-dominated state in a density range in agreement with 
that found in recent experiments on quantum-well heterostructures. 
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The rich physics arising from the Goulomb attraction 
between electrons and holes in layered semiconductor sys¬ 
tems continues to generate fundamental and technolog¬ 
ical interest. Collective many-body effects, such as the 
Fermi-edge singularities (FES) in absorption spectra pre¬ 
dicted by Mahan [l[, dominate at high carrier densities, 
while excitonic species form in dilute systems Q. Neu¬ 
tral excitons, consisting of bound electron-hole pairs, and 
charged trions, which are bound states of two electrons 
and one hole, or two holes and one electron, are elemen¬ 
tary quasiparticles that can be created via photoexcita¬ 
tion or by chemical or electrical means in a wide range 
of materials. Many optoelectronic devices, from photo- 
voltaics and light-emitting diodes Q to optoelectronic 
storage devices 0 , interconnects and switches d, @ , ex¬ 
ploit excitonic effects, as does nature in, for example, 
plant photosynthesis. 


Modulation doped or gated semiconductor quantum 
wells, including GaAs wells with AlGaAs barriers and 
InGaAs/InAlAs junctions, offer convenient experimental 
access to high-mobility electron gases. Control of carrier 
density and creation of electron-hole pairs via photoex¬ 
citation are readily achieved and early studies identified 
the FES 0, i at high carrier densities. Kheng et al. 
0] identified negatively charged trions in CdTe quan¬ 
tum wells at low carrier densities, 35 years after their 
prediction by Lamport 0- Particle confinement in a 
quasi-two-dimensional well increases the overlap of the 
hole and electron wave functions and increases the bind¬ 
ing energy of the hole compared to the bulk semiconduc¬ 
tor ll|, making this an ideal system in which to study 
trions. The crossover between high-density FES domi¬ 
nated by many-body correlations, which we refer to as 
the collective excitonic state, and the low-density behav¬ 
ior, characterized by the presence of excitons and trions, 
has been investigated experimentally 0, [12-0 . 

Huard et al. 0 observed a gradual change in the ab¬ 
sorption spectrum of modulation doped CdTe semicon¬ 


ductor quantum wells from discrete excitonic peaks at 
low carrier density to broad FES at high density. Rapid 
changes in line shapes and transition energies seen in 
the absorption and photoluminescence spectra of a gated 
modulation-doped GaAs quantum well allowed Yusa et 
al. [l2| to locate a critical “crossover” density. Simi¬ 
lar methods were used by Bar-Joseph 0. Yamaguchi 
et al. [0 recently showed that the photoluminescence 
linewidth of a gated, undoped GaAs quantum well as a 
function of energy shift in a perpendicular electric field 
can be used to measure the spatial extent of the trion 
at a given carrier density. The measured trion radius 
increases sharply above a critical density, which is iden¬ 
tified as the crossover. A theoretical description of the 
crossover was given by Hawrylak [l5| via an approximate 
treatment of electron-electron interactions. 

Despite the wealth of experimental information, inter¬ 
pretation of spectroscopic data is often not straightfor¬ 
ward and the properties of excitonic states are still de¬ 
bated. Emission of photons at the anticipated exciton 
frequencies does not unambiguously signal the presence 
of excitons, as a system formed by an electron gas and a 
hole will resonate at the exciton frequency due to many- 
body interactions [0. In addition, experimental sam¬ 
ples exhibit great sensitivity to variables such as tem¬ 
perature, finite quantum-well width, and the presence of 
disorder and localization effects from the, albeit spatially 
removed, dopants in modulation-doped quantum wells. 
Further theoretical insight into such electron-hole sys¬ 
tems is urgently needed. We have therefore performed 
variational and diffusion quantum Monte Carlo (VMC 
and DMC) 0 calculations to understand the important 
limit of a zero-temperature, two-dimensional (2D) system 
comprising a single hole immersed in a 2D homogeneous 
electron gas (HEG) interacting via the Coulomb (1/r) 
interaction. 

The relevant length scale in a 2D HEG is its density 
parameter Vg = lly/Tvn, where n is the number density. 
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while the excitonic length scale is the exciton Bohr radius 
Og = 47reoe/i^/where /i = memh/{rUe + ruh) is the 
reduced mass of the electron-hole pair, rrie and rufi are 
the electron and hole effective masses, and e is the static 
dielectric constant of the host material. The energy scale 
of excitonic systems is the exciton Rydberg (Ry*), where 
IRy* = /re^/(327r^eo£^/i^). We use Hartree atomic units 
{h = \e\ = me = dTregc = 1) unless otherwise stated. 

We use the casino code [l^ and VMC and DMC 
methods to simulate systems containing 86 electrons and 
a single hole in a periodic cell in the presence of a uniform, 
neutralizing background charge density. In VMC, expec¬ 
tation values are evaluated using a trial wave function 
containing optimizable parameters. The more accurate 
DMC method projects out the lowest energy state with 
the same nodal surface as the trial wave function 0 - 
The accuracy achieved for expectation values of opera¬ 
tors that commute with the Hamiltonian is determined 
by the trial nodal surface, while the accuracy of other 
expectation values and the statistical efficiency achieved 
are influenced by the quality of the entire trial wave func¬ 
tion. Quantum Monte Carlo (QMC) methods have previ¬ 
ously been used to study related systems including three- 
dimensional (3D) electron-hole gases l2(i. 2D electron- 
hole gases [m , trions in 2D material s l22L [ 2 ^ , excitons 
and biexcitons in bilayer systems 2^ 26l|. and positrons 
immersed in 3D electron gases 27. l28|. 

We use a trial wave function of the form 


^T(R) = e''(^)^s[X(R)] , (1) 

where R denotes the particle coordinates, is a Jas- 
trow factor that describes electron-electron and electron- 
hole correlations [ 2 ^, X(R) is a set of backflow- 
transformed coordinates [3lL l32| , and 


'l>s(R) = det 





det 




( 2 ) 


is a product of Slater determinants containing orbitals 
that pair each electron with the hole, where rj is the po¬ 
sition vector of the jth electron of spin a and is the 
position of the hole. We use a novel form of flexible pair¬ 
ing orbital whose parameters are optimized within VMC 
and which provides an accurate description of electron- 
hole and electron-electron correlation. 


(t>i{v) = exp [uGi{r)] exp{zGi • [r - r]G,ir)]r} , (3) 


where f is the unit vector in the direction of r. The 
orbital-dependent electron-hole Jastrow function uCi and 
orbital-dependent electron-hole backflow function 770 ^, 
where Gi is the ith shortest reciprocal lattice vector, 
are described in the Supplemental Material 33|. The 
free parameters in the backflow transformation, Jastrow 
factor, and orbitals are optimized using energy minimiza¬ 
tion techniques . We impose the electron-electron and 
electron-hole Kato cusp conditions 3^ via the Jastrow 
factor. 


The pairing orbitals lower the VMC energies, roughly 
halving the difference between the VMC and DMC relax¬ 
ation energies, compared to equivalent calculations with 
plane-wave orbitals. The DMC energies are only lowered 
slightly by the use of pairing orbitals; this insensitivity to 
changes in the nodal surface indicates that the effect of 
the fixed-node approximation on the energies is small. A 
detailed comparison of results obtained using optimized 
orbitals and plane-wave orbitals can be found in the Sup¬ 
plemental Material [s^. 


We use full shells of electrons in hexagonal simulation 
cells subject to periodic boundary conditions. In our pro¬ 
duction calculations we use systems with Ne = 86 elec¬ 
trons, which we deemed sufficiently large after a finite- 
size-effect analysis involving systems of up to W = 146 
electrons 0. We use a cell area of (Ne — IjTrr^, so 
that the electron density far from the hole is correct 
Mill- We calculate the electron-hole relaxation energy. 


also called the electron-hole correlation energy, by sub¬ 
tracting the energy of a HEG of the same area containing 
the same number of electrons. The energy required to 
create an electron-hole pair, for example via photoexci¬ 
tation, is given by the sum of the band gap, Fermi en¬ 
ergy and relaxation energy, the latter arising from the re¬ 
sponse of the electron gas to the point-particle impurity. 
Finite-size effects in the electron-hole relaxation energy 
are small, and the pair-correlation functions (PCFs) are 
well-converged with respect to system size . We study 
systems with mass ratios mh/me = 0.5, 1, 2, 4, and 8. 


The calculated electron-hole relaxation energies are 
shown in Fig. [T] together with the energy of the iso¬ 
lated exciton. Ex = —4 Ry* 0. We fit the electron- 
hole relaxation energies for each mass ratio to functions 
that tend to the energy of an isolated negative trion at 
Ts/ciq —^ 00 , which is the low carrier density limit of 
the electron-hole relaxation energy. We evaluate isolated 
trion energies in separate DMC calculations; numerical 
values are shown in Fig. and tabulated in the Sup¬ 
plemental Material [s^. Since trions are composed of 
inequivalent particles, the wave function is nodeless and 
DMC is exact in this case. To obtain the /i/we —>■ 0 
limit, we minimize the DMC energy as a function of the 
separation of two fixed electrons, finding an equilibrium 
separation of 0.51454(2) Og. Our data are in good agree¬ 
ment with previous results [110 in the limiting cases 
studied in those articles, and are compatible with recent 
calculations using the l/r potential for trions in transi¬ 
tion metal dichalcogenides [23| , although our energies are 
somewhat different to those of Ref. R|. 


At low densities the localized trion is similar to an elec¬ 
tron, and therefore the variation of the relaxation energy 
with Tg is dominated by the energy required to remove 
one electron from the HEG, which is minus its Fermi en¬ 
ergy. At high densities the localized trion does not form, 
and at T's/og 0 the relaxation energy diverges towards 
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FIG. 1. (Color online) Electron-hole relaxation energies for 
mass ratios mhlrrie. = 0.5, 1, 2, 4, and 8, and HEG densities 
Ts = 1, 2, 4, 6, and 10 a.u. The brown dashed line indicates 
the isolated exciton energy at —4 Ry*, and the magenta dot¬ 
ted, black short-dashed, red long-dashed, green dot-dashed, 
and blue dot-dot-dashed lines show isolated trion energies at 
mass ratios rrih/m^ = 0.5, 1, 2, 4, and 8, respectively. Fits of 
the electron-hole relaxation energies at each mass ratio have 
been constrained to tend to the respective trion energies at 
Ts/aQ 00 . Error bars are smaller than the symbols. 
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FIG. 2. (Color online) Total energy of an isolated trion as a 
function of [ 1 ,/me. (circles). The energy of the neutral exciton 
is also shown (dashed line). Error bars are smaller than the 
symbols. 


We calculate PCFs g(r) using extrapolated estimation 
[H, li^ to eliminate the leading-order dependence on the 
trial wave function. The VMC and DMC PCFs calcu¬ 
lated using plane-wave or optimized orbitals are almost 
identical, indicating the robustness of the results. Repre¬ 
sentative PCF data are shown in the Supplemental Ma¬ 
terial [13 . In Fig. [3] we show integrated electron-hole 


PCFs which give the total electron weight within a cir¬ 
cle of radius r centered on the hole. As the density is 
lowered, the behavior associated with the electron gas 
is modified by the formation of a plateau at a weight of 
two electrons, indicating the emergence of a trion. Fig¬ 
ure [3] shows that the trion radius decreases somewhat for 
heavier holes; in excitonic units (not shown) these curves 
coincide for r <rs. The electron-hole PCF tends to that 
of an isolated trion at low density, while at high density 
the free electrons screen the attractive potential from the 
hole, preventing trion formation. An intuitive explana¬ 
tion of the trends seen in the relaxation energy is afforded 
by its close relationship with the electron-hole PCF. At 
low density, the relaxation energy approaches the sum of 
the trion energy and the cost of the reduction in density 
of the surrounding electron gas far from the trion caused 
by the addition of the hole. For Vg ^ a q, the screening 
charge becomes more localized close to the hole, increas¬ 
ing the relaxation energy at high density . 



FIG. 3. (Golor online) Cumulative (integrated) electron-hole 
PCFs normalized to show the total electron weight within a 
circle of radius r centered on the hole. Data are plotted for 
constant mass ratio mh/rae = 1 and HEG densities Va = 1, 2, 
4, 6, and 10 a.u. (left panel), and for constant HEG density 
Ta = 6 a.u. and mass ratios m^/me = 0.5, 1, 2, 4, and 8 (right 
panel). 


The on-top electron-hole PCF geh{^) is proportional 
to the rate of electron-hole recombination. It could also 
be used to create semilocal two-component exchange- 
correlation functionals for use in density functional the¬ 
ory calculations for modeling holes immersed in inho¬ 
mogeneous 2D systems [i^. Considering the limits of 
an exciton in a dilute electron gas and a trion without 
electron-electron interaction in a dilute electron gas, we 
propose a relation [s^: 

5e/i(0) = c/i^rf-b 1 , (4) 

where c is a dimensionless parameter that is roughly in¬ 
dependent of [1 and Tg and takes values slightly above the 
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exciton limit of c = 8. We have extrapolated the PCFs 
to r = 0 and plotted the results against fiVg (= Tg/aQ) 
in Fig. 01 Equation (0]) fits the data well over the pa¬ 
rameter space studied. We obtain c = 9.742(7) from the 
fit; the variation of c with /r and rg is analyzed in the 
Supplemental Material (33| . 

The value of the electron-hole PCF at r is the ratio of 
the electronic density a distance r from the hole to that of 
the surrounding electron gas. Thus, the value of the PCF 
at its first minimum, also shown in Fig. 01 measures the 
degree of isolation of the localized trion. The minimum 
in the electron-hole PCF develops at about Vg ~ Og and 
rapidly becomes more pronounced as the carrier density 
decreases, with the minimum electron density falling be¬ 
low 25% of that of the surrounding HEG by rg ~ 4.5 Ug. 




FIG. 5. (Color online) Crossover of hole-in-HEG system 
from the high-density collective excitonic state to a localized 
trion immersed in a low density HEG as a function of mass ra¬ 
tio mh/rrie. and density parameter rs/ag. Experimental data 
showing the evolution of absorption and photoluminescence 
spectra from the Fermi-edge singularity to discrete trion and 
exciton peaks are shown as colored areas, and are consistent 
with our results. 


mize noise during optimization. In addition, each result 
is an average over eight independently optimized wave 
functions [33[. We have applied this method to compute 
the momentum densities of four representative systems, 
as shown in Fig. 01 


FIG. 4. (Color online) Values of the electron-hole PCF at 
r = 0 (on-top PCF) and at its first minimum. The solid line 
is a least-squares fit of the on-top PCF data to Eq. (0]). The 
PCF values 0.25, 0.5, and 0.75 are represented with dotted 
grey lines. Error bars are smaller than the symbols. Note the 
logarithmic scales. 


The gradual emergence of trions in our QMC calcu¬ 
lations is compared to experimental data from semicon¬ 
ductor quantum-well systems in Fig. [5l This gradual 
crossover occurs in a parameter range consistent with 
the absorption and photoluminescence spectra seen ex¬ 
perimentally. The values of the dielectric constant and 
effective masses used are given in the Supplemental Mate¬ 
rial, together with an alternative picture of the crossover 
in the same parameter range [s^. 

We have calculated the electron-hole center-of-mass 
momentum density p{k) by constraining one electron to 


remain on top of the hole [2^ [SSj, [STJ. Since it is not 


possible to use extrapolated estimation for this quantity, 
we report VMC results using the optimized orbitals. The 
momentum density is sensitive to the quality of the trial 
wave function, which has been compensated for by us¬ 
ing 500,000 configurations, a much larger number than is 
required to converge the VMC energy, in order to mini- 



FIG. 6. (Color online) Momentum density for four representa¬ 
tive points of our “phase diagram”: left column, rtihlme = 1; 
right column, mhlrrie, = 8; top row, = 1 a.u.; bottom 
row, Ts = 6 a.u. Error bars are shown but are sometimes 
smaller than the symbols. Our results are normalized such 
that 2nkp(k) dk = nk%. 


At high density we obtain strong enhancement of mo¬ 
mentum density just below the Fermi edge, together with 
a small tail above the edge. This behavior was predicted 
theoretically by Carbotte and Kahana [4^ 14^ and re- 
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cently demonstrated numerically by Drummond et al. 
[28| for a positron in a 3D HEG at metallic densities. 
However, our results for 2D systems also show the for¬ 
mation of a small, broad peak above the Fermi edge as 
the density is lowered. In contrast to the momentum 
density for k < kp, this unusual peak is insensitive to 
small changes in the wave function parameters and to 
the precise form of optimizable wave function used. The 
peak emerges gradually as the density is lowered from 
that at which we estimate trion formation to begin, and 
becomes higher and narrower, its center moving closer 
to the Fermi edge, as the density is lowered to rg = 10 
a.u. We demonstrate that this peak is associated with 
the formation of a trion in the Supplemental Material 
0 . 

In conclusion, we have performed highly accurate 
QMC calculations for a system containing a single hole 
immersed in a 2D electron gas. Our results demonstrate a 
crossover between a collective excitonic state and a trion 
state as the density of the electron gas is lowered and as 
the mass ratio is increased. The electron-hole relaxation 
energy, PCF, and electron-hole center-of-mass momen¬ 
tum density each show evidence of the crossover. The 
density and mass range in which trion formation begins 
is in good agreement with recent experiments [1, 12, l3 ]. 
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